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Abstract 

We consider the Navier-Stokes equation on EI^(— a^), the two dimensional hyperbohc space 
with constant sectional curvature —a^. We prove an ill-posedness result in the sense that the 
uniqueness of the Leray-Hopf weak solutions to the Navier-Stokes equation breaks down on 
BI^(— a^). We also obtain a corresponding result on a more general negatively curved manifold 
for a modified geometric version of the Navier-Stokes equation. Finally, as a corollary we 
also show a lack of the Liouville theorem in the hyperbolic setting both in two and three 
dimensions. 

1 Introduction 

We investigate the impact the geometry of the underlying space has on the Leray-Hopf solutions to 
the Navier-Stokes equation. More precisely, we consider the Navier-Stokes equation on negatively 
curved manifolds and present how the negative scalar curvature causes the break down of the 
uniqueness of the Leray-Hopf solutions. 

Before we state the main results, we survey some necessary historical background from both 
geometric analysis, and the regularity theory for the Navier-Stokes equation in the usual Euclidean 
setting. 

1.1 Regularity theory for the Navier-Stokes equation on M" 

The Navier-Stokes equation on the Euclidean space R" is given by 



dtu- Au + u-Vu + VP = 0, 
divu = 0. 



(N-Sm.) 



Long time ago, for the dimensions n = 2 and rt = 3, Leray FblJ and Hopf [7] established the 
existence of global weak solutions u £ L°°{0, oo; L2(]R")) n 2.^(0, oo; iji(E")). Due to their work, 
we now have the following general existence result, which historically served as the foundation for 
further works in the regularity theory for N-Sr^. 

Theorem (Leray-Hopf weak solutions [HI [7]). Given any initial datum uq E iy^(K") . there exists 
at least one -valued function u G L°°{0,oo; L'^{W)) H L^{0,oo; (&.""■)) which satisfies the 
following properties 
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• For any (j) = (0i, ...,(/)„) G C^((0, oo) x M") with div (f> = 0, we have 



f 

Jo 



■ oo 



■ dt(p + '^{d.j4>i){djUi) - y^^{dj(pi){uiUj)dxdt = 0. 




• For every t > 0, u satisfies the following global energy inequality 




\u{t,x)\'^dx + 2 / iVupdxds < / \uo\'^dx 




• it(0, •) coincides with the initial datum uq in the sense that 



limt^Q+ \\u{t, •) - -uo||l2(R") = 0. 



(1.1) 



Now that we have existence of the Leray-Hopf solutions for N-Sr2 and N-Sr3, we proceed to 
address the question of regularity. The regularity of Leray-Hopf solutions on greatly differs from 
the corresponding regularity problem for Navier-Stokes equation on R^. Indeed, the smoothness 
and uniqueness of Leray-Hopf solutions for N-Sh2 was established in the work of Leray (see for 
instance f22J). As a sharp contrast, the regularity and uniqueness of solutions to the N-Sr3 equation 
is a long standing open problem although due to the concentrated efforts by generations of PDE 
specialists there has been a significant progress in this area. 

Because of the limitation of space, we only mention some typical regularity criteria for Leray- 
Hopf solutions to N-Sr3. We also note that one of the goals of this discussion is to illustrate why 
there is more focus on the question of regularity than that of the uniqueness. 

Now, the first significant effort to break the silence since the fundamental work of Leray and 
Hopf, was made in 1960's through the efforts of Prodi [17], Serrin [20], and Ladyzhenskaya [T2] 
leading to the following regularity and uniqueness result (for more historical remarks see for 
instance [5]). 

Theorem 1.1. [Prodi, Serrin, Ladyzhenskaya] Let u e L°°{0,T; L'^{R^)) D L'^{0,T; H\E.^)) 
to be a Leray-Hopf weak solution to N-Sr3, which satisfies the additional condition that u G 
LP{0, oo; L'^(R^)), for some p, q satisfying | -t- | = 1, with q > 3. Then, u is smooth on (0, T] x 
and u is uniquely determined in the following sense 

• suppose V e L°°(0, T; L^(M^)) n L^(0, T; 7?^(R^)) is another Leray-Hopf weak solution such 
that u(0, ■) = f (0, •). Then, it follows that u = v on (0,T] x R^. 

Here, let us briefly mention why the case of g = 3 was not included in Theorem ll.il Indeed, 
it is well known that a solution 9 : [0, T) x R"^ R to the heat equation arising from any initial 
datum 9q G L'^(R'^) satisfies the following estimate for any pair of indices p, q with | + | — 1 ^''^d 
q > 3 (see [9] Appendix]). 



where Cq depends only on q. So, in some sense, the extra condition as imposed on the Leray-Hopf 
solution u in Theorem 11.11 ensures that the qualitative behavior of the Leray-Hopf solution u would 
be a slight perturbation of solutions of the heat equation. Another explanation for the relatively 
simple status of Theorem 1 1.1 1 is that the L^L'^ norm of the solution u under the integral condition 
as promised in Theorem 11.11 shrinks to zero under the natural scaling Ug(t,x) = eu{e^t,ex) as 
e — >■ 0. However, this is no longer valid in the borderline case of i°°(0, T; L^{M.'^)). This partially 
explains the long delay in the settlement of this exceptional case of m G L°° {Q,T; L^{M.'^)), which 



II^IIlp(0,T;L9(R3)) < C'(9)II"0|1l3(r3), 



(1.2) 
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was finally established in the recent work of Escauriaza, Seregin, and Sverak [S] . 

Before we close our discussion let us mention that one of the working principles in the reg- 
ularity theory of Navier-Stokes equations on M'^ is (more or less) to reduce the situation under 
consideration (say u G L°°{L^) in the case of [S]) to the regime which is covered by Theorem ll.il 
Once this can be achieved, then the uniqueness of the solution would come for free, due to the 
uniqueness claim in Theorem 11.11 This explains to some extent the fact that regularity issue is 
more of a focus than the uniqueness issue in the regularity theory for Navier-Stokes equations in 
the M'^ setting. However, as is well-known, the weak formulation for Leray-Hopf weak solutions 
to the Navier-Stokes equation on M.^ only gives the natural bound u G ^^(0, oo; L''(R^)), with 
indices p, q satisfying - _|_ 3 _ | _ readily sees that there is a significant gap between the 

natural bound oflEcred f)y the weak formulation and the condition required by Theorem II. 1[ and 
it is unclear how to make a link between them. See again the introduction of [9] for a discussion 
about this point, and for further developments, we refer our readers to a piece of recent work by 
Vasseur [21]. 

1.2 Navier Stokes equation on a Riemannian manifold 

Historically speaking, the correct form of the Navier-Stokes equations in the Riemannian manifold 
setting was first obtained by Ebin and Marsden [4 . They considered compact Riemannian, ori- 
ented, n-dimensional manifolds both with and without boundary. Moreover, they remark that the 
derivation of the correct equations assumes that the manifold in question is Einstein, i.e.,Ric = Xg, 
for some constant A where Ric is the Ricci tensor and g is the Riemannian metric g. We note, this 
is in particular true of space forms, where Ric — {n ~ l)KMg (see Section [2.1! below) . 

According to [4] the ordinary Laplacian should be replaced by the following operator in the 
formulation of the Navier-Stokes equation on a Riemannian manifold 

L = 2Def *Def = V* V + dd* - Ric = {dd* + d*d) + dd* - 2 Ric, (1.3) 

where Def and Def* are the stress tensor and its adjoint respectively, V stands for the induced 
Levi-Civita connection on the cotangent bundle T*M, V V stands for the Bochner Laplacian, 
with V to be the adjoint operator associated to V, {dd* + d*d) = —A stands for the Hodge 
Laplacian with d* to be the formal adjoint of the exterior differential operator d, and Ric is the 
Ricci operator (see Sections 12.11 and 12.61 for definitions and [31 [22] for a further discussion of the 
deformation tensor). We first remark that the operator L as given in expression (|1.3p is an operator 
sending sections of T*M into sections of T* M . This means that, the Navier-Stokes equations on 
a Riemannian manifold M is formulated in terms of sections of T* M instead of vector fields on 
M. 

As a result, the usual convection term V^u in terms of vector fields also has to be rewritten. 
There is a natural correspondence between vector fields and 1-forms (see Section 2.1), which 
produces the term S/uU* ^ where U is the unique vector field corresponding to a 1-form U* . 

In summary, we regard the solutions to the Navier-Stokes equations on a general Riemannian 
manifold M to be differential 1-forms U* £ C°° {M;T* M) satisfying the following differential 
equation 

dtU* +L{U*) + VuU* +dP = 0, (1.4) 
d*U*=0. (1.5) 

where P is a scalar function on AI. Using the fact that U* is divergence free we can further rewrite 
the equations as follows 

dtU* -AU* +^uU* -2Ric{U*)+dP = 0, 
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which is the main equation that we study in this article. 

A less natural equation to consider is one without the Ricci operator. We refer to it as the 
modified Navier-Stokcs equations on M and record it here 

dtU* - AU* +yuU* + dP = 0. , ^ 

(1.6) 

d*U* =0 ' 

It is less natural from the point of view of the derivation of the Navier-Stokes equations performed 
in [J. We consider it in this paper, because we would like to present how a more general manifold 
than just a space form can influence the behavior of solutions (we explain this more below). 
Since Def * Def U* plays now the role of the dissipation, the global energy inequality becomes 

\U*\^{t,x)+2 [ [ -g{V)eiU*,V)eiU*)ds< I |uoP, (1.7) 

M Jo J M Jm 

where g(-,-) stands for the inner product structure on the bundle T* M (g) T*M induced by the 
Riemannian metric g{-,-) on M (see Section l2.ip . 

We now mention some of the previous results on a Riemannian manifold (for more see [3] 
and references therein). Priebe [16] appears to be the first one to follow [4] and use the correct 
version of the equations N-Sm instead of ()1.6p . [I6] also assumes compactness of M and works on 
manifolds with boundary. Dindos and Mitrea [3] consider the linearized version of the stationary 
Navier-Stokes equations on a subdomain of a compact Riemannian manifold. In fact, we have not 
been able to find any results for non-compact manifolds except for the result of Q.S. Zhang [55]. 
In !26| the author shows the ill-posedness of the weak solution with finite norm on a connected 
sum of two copies of M.^. Hence the topology of the underlying manifold seems to play a role. 
In this paper, we take a geometric point of view and also consider the dissipation term, which 
involves careful computations. 

We are now ready to state our main results. 



1.3 Statements of the results and discussion of the proofs 

Theorem 1.2 (Non-uniqueness of N-SH'2(-a2)). Let a > 0. Then there exist non-unique Leray- 
Hopf solutions to N-SH2(-a2). 

Remark 1.3. The consequence of Theorem 1 1.21 is that unlike in the Euclidean setting, the notion 
of the Leray-Hopf solutions might not be the proper foundation for the study of solutions to the 
Navier-Stokes equations on the space form with negative sectional curvature —a^ in dimension 
two. The question of what happens to the strength of the framework of the Leray-Hopf solutions 
on H^{—a?) is open. The lack of a definite answer at this point is mainly due to the specific 
form of the techniques we use in this paper (see remarks below). It is an interesting question to 
sec if perhaps the techniques could eventually be extended/modified to give some insight into the 
corresponding question in 3 dimensions. 



Corollary 1.4 (Lack of the Liouville theorem for space forms). Let n > 2, and a > then there 
exist nontrivial bounded solutions o/ N-S/fii(_Q2) . 

Remark 1.5. The proof of CoroUarv II .41 and Corollarv 11.71 below follows trivially from the proofs 
of their theorems. Moreover, it does not require any of the delicate estimates developed in this 
paper. As such it is just a by-product of the main results and we only include it here for com- 
pleteness, and because of the general importance the Liouville theorems play in the subject of the 
Navier-Stokes equations. See Section [6.31 for motivation and some background. 
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If one decides to omit the Ricci term from the equation, we can also have a non-uniqueness 
resuh on a more general negatively curved Riemannian manifold than just ]HI^(— a^). 

Theorem 1.6. Let a,b > be such that < a < b, and let M be a simply connected, complete 
2-dimensional Riemannian manifold with sectional curvature satisfying —b^ < Km < . Then 
there exists non-unique Leray-Hopf solutions to the modified Navier-Stokes equation ()1.6|) . 

Corollary 1.7 (Lack of the Liouville theorem in the hyperbolic setting). Let n > 2, and b > a > 

and let M be a simply connected, complete n-dimensional Riemannian manifold with sectional 
curvature satisfying —b^ < Km < — a^. Then there exist nontrivial bounded solutions of (11.61) . 

Remark 1.8. Note, the lower bound ifo < a is no longer required in the corollary. See the 
discussion below for why the lower bound is present in Theorem II. 6[ which also explains why we 
do not need it in Corollarv ll.7l 

The above results are based on the abundance of the nontrival bounded harmonic functions 
in the hyperbolic setting. Such abundance is ensured by the works of Anderson [1] and Sullivan 
[2T| . Our idea of trying to benefit from them was inspired by a remark of Tsai [23l Remark 
5.4]. [53] eliminates a possibility of self-similar solutions to N-Sr3 (which merely would satisfy 
the local energy inequality) by showing that existence of the self-similar solutions is equivalent 
to solving a certain stationary system. Without assuming enough decay, one could construct 
nontrival solutions of the system in question in the form U — VF, and P = —^\U\'^ — ay ■ U, 
where F is a harmonic function on M.^ and a > 0. In our case, due to [T] and [2T] we have a 
plethora of nontrival bounded harmonic functions, which gives us a basis for this article. 

The solution pairs {U* , P) we consider have the following form 

U* = ^{t)dF, 

P = -dtTpit)F -^\dF\^ + 2a^F, ^^'^"^ 

where ipit) = exp(— ^) for some A > 2c? , and is a nontrival bounded harmonic function on 
H^(— a^). Verifying that {U* ,P) solves N-Se2(_o2) is simple when we use Hodge theory (see Sec- 
tions and [S]) and Lemma [Ol 

Before we proceed any further, we remark here that the differential geometric work in [Tl 1211 12] 
ensures the existence of nontrival bounded harmonic functions on a more general negatively curved 
Riemannian manifold with suitable lower and upper bounds imposed on the sectional curvature. 
On the other hand, the existence of nontrival bounded harmonic function on H"(— 1) is an old 
classical result obtained through an integral representation formula with an explicit Poisson ker- 
nel on the Poincare ball model for the space form ]EII"( — 1) (for more details see the work of Hua 
[8]). However, such classical approach relies heavily on the explicit formula of the Poisson kernel 
derived from the group of isometrics of the space form !!"(— 1). It seems that, as compared with 
the differential geometric approach of [T] dT] [2] , such classical approach does not reveal the role 
played by the negative sectional curvature of the hyperbolic manifold in the existence of nontrival 
bounded harmonic functions on H"(— 1). 

The last remark may explain why the proper generalization of the above mentioned classical 
result to the more general setting of negatively curved Riemannian manifolds was only established 
in the more recent works [TJ [21] [2]. Since we intend to show not only the break down of the 
uniqueness of Leray-Hopf solutions in the hyperbolic space setting, but more importantly the de- 
cisive role played by the negative sectional curvature of a hyperbolic manifold in causing such a 

^In fact, taking solutions of the form ip{t)\7F seems to be a well known convention, and we just happened to 
learn about it from | 23| . But also see |26| . 
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breakdown, we will unconditionally choose the differential geometric framework as established in 
[U [211 [2] as the basic ground in this paper. 

Moreover, since the differential geometric machinery as demonstrated in [T], [3T], and [5] is 
designed to establish the existence of bounded nonconstant harmonic function on a general neg- 
atively curved Riemannian manifold which lacks the homogeneity property enjoyed by the space 
form ]HI^(— a^), the only best way to justify the use of such differential geometric machinery in our 
paper is to cast our theorems, lemmas, propositions in the most general setting of a negatively 
curved Riemannian manifold at the starting point of the paper. However, we slowly narrow down 
our setting by imposing further restrictions on our results whenever such restrictions are needed in 
proving the finite integral of a certain function or in handling the extra Ric term in the formulation 
of the Navier-Stokes equations. 

As stated in Theorem II .Gi our non- uniqueness result also holds for a more general negatively 
curved Riemannian manifold with the lower bound —b^ and the upper bound —a^ of the sectional 
curvature satisfying < | < a < 6, provided if the extra Ric term in the Navier-Stokes equation 
is dropped. Indeed, the final restriction to the space form ]HI^(— a^) is required only because of the 
presence of —2 Ric in the formulation of the Navier-Stokes equation N-Sa/. 

Now, we explain our strategy in establishing the finite energy and the finite dissipation of the 
time dependent velocity field U* = ip{t)dF as given in equation (|1.8p . We start our discussion 
by saying that our exposition is based on the material in the second chapter of the book |19j by 
R. Schoen and S.-T. Yau. In the first section of the second chapter of [19], one sees that, with 
the prescribed function (f> g C^{S{oo)) given on the geometric boundary S'(oo) (see Section [2^ 
attached to the n-dim complete, simply connected Riemannian manifold M with sectional curva- 
ture satisfying —b^ < Km < ~a? < 0, the bounded harmonic function F on M which satisfies 
the Dirichlet boundary condition F|5(oo) = <?!> is sought after by means of creating two barrier 
functions (j) — ae"^^ and 4> + ae~^P, which serve as the lower bound and upper bound for F and 
where p stands for the distance function on M from a selected base point O in M (also see Section 
I2.5l and l2.2|) . This is done in the spirit of the classical Perron's method. But such an application 
requires the subharmonicity of (j) — ae~^'' and the super harmonicity of 0-1- ae~^'' , whose validity 
critically depends on the following two facts (for details, see [T^ ) 

• Laplace comparison theorem: If Km < ^a^, then, the Laplacian of the distance function p 
(from a selected base point O in M) satisfies Ap > (n — l)acoth(ap) > (n — l)a. 

• the smooth function (j) is constructed in a specific way so that we have (j)\s{oo) — 4' and that 
the oscillation of over any geodesic ball Bx{l) in M has exponential decay of order e^°'''^^\ 
for any x £ M . 

Due to the above two facts, it can be deduced that the choice of the (5 > 0, which ensures 
A[(/) — ae"^''] > and A[(j) + ae~^P] < (and hence the success of the Perron's method), has to 
satisfy the constraint S < a (see Section r2. 5. ip . 

Based on what we learn from the above construction of the bounded nontrival harmonic func- 
tion on M we employ, in Section 3 the gradient estimate for harmonic functions due to S.-T. 
Yau [25] to show that the decay rate for |VF|(x) as p{x) approaches infinity is at least of the order 
Q-Spix) ^ for any S < a. That is we have 

\\/F\<C{a,S)W\\sia.)e-'', 

on M. Here, we want to mention that, with the hindsight from one of the two Harnack's inequalities 
as established in the second chapter of jlj] , one can argue that such an exponential decay for the 
gradient of our bounded harmonic function is more or less expected and may not be surprising. We 
believe that such an exponential decay could be more or less well known to researchers working 
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in geometric analysis. But in any case, we give a clean and simple proof of it in Section 3. 

Next we note that the exponential decay |Vi^| < C{a,d)\\(j)'\\s{oci)P'^^'' not only gives, in the 
special case of the two dimensional space form ]HI^(— a^), the -finite property of | VF| on H^(— a^) 
(and hence the finite energy property of the velocity field U* = ■tp{t)dF as given in (ll.Sp ). but also 
demonstrates the limitation which prevents us to draw the same -finite property of \'S/F\ in the 
setting of the three-dimensional space form H'^(— a^) . This limitation mainly comes from the fact 
that oscb^(i) (j) only has exponential decay of order e~°'^'^^\ which prevents us from choosing a 
6 > larger than a in the Perron's method; yet the growth rate of the volume form on ]HI^(— a^) is 
exactly of the order ^ sinh^(ap). For the first time, we encounter an obstacle which forces us to 
restrict our theory only to the case of 2 dimensional negatively curved Riemannian manifold M 
with < j^^j < _q2 ^ 

We observe that, up to this point, the lower bound condition Km > —b^ has not been involved 
in the big picture yet. However one does eventually have to pay a special attention to the relative 
largeness of b when compared with a since the lower bound Km > —b^ of the sectional curvature 
determines the growth rate of the volume form of M via the comparison theorem for Jacobi fields 
in differential geometry. More precisely, with the condition Km > —b^ imposed, the growth rate of 
the volume form of the 2-dim negatively curved manifold M is at most ^ smh{bp). Yet, again, the 
decay rate of |VJ^| is of the order e"*", with any S < a. As a result, the survival of the property 
\\^P\\l^{m) < oo critically depends on the competition between the decay rate e~^°' of |Vi^| and 
the (possible) worst growth rate sinh(5p) of the volume form of M. This fully explains the need 
to further restrict our setting by imposing the condition < | < a < &, so that the parameter 
S can fit within the range of ^ < S < a, which is enough to ensure the survival of the -finite 
property of \S7F\ on the 2-dim negatively curved manifold M. 

Once the L^-finite property of \\7F\ is established for 2-dim Riemannian manifold satisfying 
—b^ < Km < —a^ and < | < a < 6, we proceed to show the finite dissipation of U* = 4'{t)dF 
under the same setting in Sections 4 and 5 of our paper, which are the most delicate parts of our 
work. In regard to this, our basic idea lies in the structure of the following important formula in 
differential geometry [19] 

A(|Vi^P)(a;) = 2g(y{VF),V{VF)){x)+2'^d^F{x)d^{/\F){x)+2Ric{VF,VF){x). (1.9) 

The formula is obtained by performing a calculation with respect to the normal geodesic coordi- 
nates about the selected point x in our 2-dim negatively curved Riemannian manifold M . Since 
A(|VFp) = div{V{\VF\'^)}, if we can show that the L^-norm of |V|VFp| is finite, then, we 
will immediately have /j^A(|VFp) = 0. Hence, it follows from the above formula that (see 
Proposition 5.1 of Section 5 for the technical details) 

/ 5(V(VF), V(Vi^)) = - / Ric{VF,\/F) <b^ I |Vi^p, (1.10) 

which gives the finiteness of the dissipation term Jj^^g{VU,VU) = [^'(t)]^ /^.^ 5(V(Vi^), V(Vi^))) 
as required in the Leray-Hopf formulation. 

Next, the required L^-finite property of | V| VFp | is established with the assistance of a covering 
Lemma 4.1. Due to the fact that Ric(VF,Vi^) > -b^jVi^p and that AF = 0, formula pTg)) 
ensures the subharmonicity of |Vi^p -I- Ae~'^^P on {x e M : p{x) > (with ^ < S < a), for 

some sufficiently large A > and R{S) > 0, both dependent on 6, which in turns allows us to 
obtain (see the proof of Proposition 4.4 for the technical details) 

• the integral of |V{|Vi^|2 + Ae-^^P}\ over any geodesic ball i?x(3(l + -j)) lying within the 
outer region {x £ M : p{x) > -^^j^} is bounded above by C(a, &; ||0||oo) e~^^'''-^\ 
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Next, since | V{| Vi^|^ + '^^p}\ is continuous on M, to see the extent to which the above fact can 
ensure the finiteness of the integral |V{|VFp + yle^^*''}|, we just have to further decompose 
the outer region {x £ M : p{x) > -^^j^} into a countable number of rings {x £ M : k — 1 < p{x) < 
k + 1}, indexed by sufficiently large positive integers k; and argue, as in our covering Lemma 4.1 
that, due to the lower bound —b^ on the sectional curvature Km, it is sufficient to use a total of 
[ttc'''^] + 1 geodesic balls with radius 3(1 + ■^) to cover the ring {x £ M : k ~ 1 < p{x) < fc + 1}, 
which in turns ensures that (see the proof of proposition 4.4 for more details) 

• the inequality /^^(g^^^i^^ |V{|VF|2 + Ae-^^P}\ < C(a, 6; ||<?!)||oo) e^^M^^) for any geodesic 

ball Bx{3{l + ^)) C {x £ M : p{x) > ^^^} is strong enough to deduce that the integral of 
|V{|VFp + Ae^'^^P}\ over the outer region {x e M : p{x) > is finite (thanks to the 

condition < | < a, which allows S to be within the range | < S < a). 

The above observation gives /^^ |V{|VF| + Ae'^^P}\ < oo. Due to the fact that J^e^'^^P < 
oo, which is ensured by the condition 2S > b (see inequality ()4.18p ). we finally conclude that 
Im ^ which is a backbone ensuring the correctness of equation (jl.lOp . 



1.4 Organization of the article 

In order to make the paper as self-contained as possible in Section [5] we collect some facts from 
the differential geometry and in particular some background specific to the negatively curved 
manifolds. 

Once that is done, we are ready to establish fundamental statements needed for the proof of 
Theorems [121 [ESI and their corollaries. They are: 

1) Exponential decay of the gradient of bounded harmonic functions on negatively curved 
manifold-Section 3. 

2) Finiteness of ||V |Vi^|^ ||ii(]a2(_a2))-Section 4. 

3) Global energy inequality tools-Section 5. 
Section 6 contains the proofs of the main results. 
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2 Preliminaries 

In this section we gather all the necessary tools from the literature needed in our proof. A lot of 
it relies on [TH], and we list it here for the convenience of the reader. We also develop some precise 
statements regarding the volume forms on the negatively curved manifolds. 
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2.1 The Levi-Civita connection, deformation tensor and Ricci curvature 

Here we recall sorae general background from Riemannian geometry (see for example [10) . |15] . 
[T5]). In particular, we take a closer look at the deformation tensor mentioned in the introduction. 

Let M be an n-dimensional complete Riemannian manifold, and TM and T*M tangent and 
cotangent bundles on M respectively. Let g be a Riemannian metric, g{-^ •) G C°°(M, S'^T*M), 
where S'^T*M denotes symmetric bilinear forms on TM, and V the Levi-Civita connection on 

V : C°°{M,TM) C°°{M,T*M (g)TM). 

Let (•, ■)tm®T' m be the natural pairing between vector fields and 1-forms on Ad . Given a vector 
field X e C°°{M- TM), using the metric g we can define X* G C°°(M; T* M) by 

{Y,X*)tm^T'M = g{X,Y), Y e C^{M-TM). (2.1) 

Similarly, given a 1-form uj £ C^{M;T*M) we can define S C^{M;TM) by 

9{Vu,Y) = {Y,L0)TM(g,T-M- 

Therefore, the Riemannian metric g gives rise to the natural identification C°°{M; TM) — C°°{M; T*M). 
In particular, if is a smooth function on M, and d is the exterior derivative we have 

(S/F)* = dF. (2.2) 

Next g also induces its dual metric g*{-, ■) G C°°{M, S'^TM) by 

9*{^,l) ^ 9{vlj,Vj), w,7eT*Af. 

Then note 

\dF\^ =g*idF,dF)=g{VF,VF) = \VFf . (2.3) 

Now by using g again, we can also induce the corresponding positive definite inner product 5(-, •) 
on the bundle T* M T* M, which is precisely characterized by the following condition 

• Let ei, 62, e„ to be a local orthonormal moving frame for TM, and let 9^, 9^, ...0" to be 
the corresponding dual frame for T*M , then, the list {6'' ^ 9^ : 1 < j,k < n} is orthonormal 
with respect to g(-, •). 

Next, the Levi-Civita connection V on the tangent bundle TM induces the associated Levi-Civita 
connection V on T* Ad by 

V:C°°iM,T*M)^C°°(M,T*M(g)T*M), (2.4) 
(Vxuj){Y) = {VxVu.)*{Y), oj e T*Ad,X,Y € TM. (2.5) 

Notice, for the sake of convenience and keeping with the conventions, we use the same abbreviation 
V to denote both the Levi-Civita connection on TM and the induced connection on T*M. No 
possible confusion should arise since the meaning of the symbol V will be clear from the context. 
In particular, we have, by the definition of the induced connection V on T* Ad , the property that 

VxY* = i^xY)*, (2.6) 

for any smooth vector fields X, Y on Ad. 

We now turn our attention to the deformation tensor 

Def : C°°{M,TAd) ->■ C°°{Ad, S^T*M), 
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defined by (see for example [3]) 

{DefX){Y,Z) = ^{g{VYX,Z)+g{WzX,Y)}, Y, Z e {M,TM). (2.7) 

Using the natural identification of the space of vector fields with the space of 1-forms on M dis- 
cussed above, the operator Def can be regarded as the operator from C°° {M, T*M) to C°° (M, S^T*M), 
and can be defined alternatively as follows: 

Definition 2.1. For any I- form uj G C°°{M,T*M), the deformation tensor Dei 6 G {M, S'^T* M) 
is given by 

(Def w)(r, Z) = ^{{Z,VY^)TAmT'M + {y,^Z^)TAmT'M}, (2.8) 

for any Y, Z e C°° {M, TM) . 

In the sequel we also need the following. If we express uj G C°°{M,T*M) locally as uj = 
ujadx"-, then Def w can locally be expressed as 

Defuj = ^ '^{'^j;k + uJk-j)dx^ (8) dx'', (2.9) 

where ujj-k = dkUJj — r'feWi, where F^^ are the Christoffel symbols. 

Before we go to the next subsection, we briefly discuss the Ricci curvature on a complete n- 
dimensional Riemannian manifold M. Recall, the Ricci curvature is a symmetric tensor Ric S 
C°°{M,S^T*M) defined as follows 

RiCp(X,X)= KM{X,e,), peM, 

i<i<ji-i 

where ei, 62, ... e„_i are some unit vectors in TpM such that {X, ei, ...e„_i} forms an orthonormal 
basis for TpM. In many occasions, we use the symbol Ric(M) for Ric. Moreover, if we write 
Ric(M) > -b^, it means that Ric(X,X) > -b^\X\^, for all X G C'^{M,TM). Moreover, it is 
clear that, for 2-dimensional Riemannian manifold M, the notion of Ricci curvature Ric coincides 
with the sectional curvature Km- 

Besides the Ricci curvature tensor Ric G C°° {M,T* M eg) T*M), we also need to consider the 
Ricci operator Ric : C°°(M, T*M) -J> C°°{M, T*M) sending the space of 1-forms into itself, which 
is defined by 

Ric(M*) = ^ ■q''g{R{ea, eb){eb), u), 

a,b 

where R{-, •) is the Riemannian curvature tensor, ei, 62, ...e„ is a local orthonormal moving frame 
for TM, and 77^, 77^, ...77" stand for the associated dual frame for T*M with respect to ei, 62, ...e„. 

In the case of the space form IHI"(— a^) with sectional curvature — a^, we have R{ea, eb){eb) = 
— Ca, for any local orthonormal frame ei, 62, ...Cn 

of TH"(— a^). Hence, in particular we have 

the following fact 

Ric(u*) = 77"g(i?(ea, eb){eb), u) 

a,b 

= -a^{n-l)J2v''9{ea,u) 

a 
a 

= ~a^in-l)u*. 
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We end this section with a quick summary of basic facts about the Ricci curvature Ric G 
C°°{M,T*M(g)T*M) and the Ricci operator Ric : C°°{M,T*M) C°°{M,T*M). In particular, 

RiCp(X, X) = Km{p) , peM,X e C°°{M, TM), dimM = 2, (2.10) 
Ric(a;) = -(n - 1)0^0;, u; e iH''{-a^),T* H'\-a^)),a > 0,n > 2. (2.11) 

2.2 Estimates and identities used 

As usual, we start with a complete n-dimensional Riemannian manifold M, and consider the 
geodesic normal coordinates on M about a selected base point O. One of the fundamental prop- 
erties of the normal coordinates, which we use in computations, is that the Christoffel symbols all 
vanish at O: 

T)k - 0, (2.12) 

(see for example |13) for more on normal coordinates). 

In the case of a complete, simply connected, n-dimensional Riemannian manifold M with 
sectional curvature satisfying —b"^ < Km < < 0, the Cartan-Hadamard theorem ensures that 
the geodesic normal coordinates on M about any selected base point O G M must be globally 
defined, which also implies in particular that M is diffeomorphic to R". Moreover, in this case, 
between any two points p, q in such a Riemannian manifold M, the geodesic joining p and q is 
unique, and hence we just use the symbol pq to denote the unique geodesic joining p and q, and 
\pq\ stands for the length of the geodesic joining p and q. 

Define the distance function from a point p G M to a point x by 

Pp{x) = \px\ . 

We usually omit the subscript p and simply write p{x) since the base point is clear from the 
context. 

Lemma 2.2 (Properties of the distance function [19, Ch. 1]). The distance function p{x) defined 
as above is smootr^ on all of M, where M is any Riemannian n-manifold such that the exponential 



map defines the diffeomorphism between M and R". In addition, 

|Vp|' = 1, (2.13) 

A/9 = (n — l)fccothfcp if M has constant sectional curvature = —k"^, (2-14) 

A/9 > (n - l)a • coth(a/9) > (n - l)a if Ric(M) < -a^ (2.15) 

A/9 < (n- 1)6- coth(6/9) < (n- l)ii^ i/ Ric(Af) > -(n - 1)6^ (2.16) 

P 



Lemma 2.3. p. 35] Let M be a n-dimensional simply connected, complete Riemannian mani- 
fold with sectional curvature satisfying —b^ < Km < — a^. Let O £ M to be a selected based point, 
and let xi, X2 be two points in M for which \Oxi\ = \Ox2\ — R, for some R > 1. Moreover, let 
9 = A{Oxi, 0x2)- Then there exists a sufficiently large universal constant Rq > 1, depending only 
on a, b, and n such that, whenever \Oxi \ = \Ox2 \ — R> Rq, we have the following 

2R+-{log9-l) < \xTx^\ < 2R+^{\og9 + l). (2.17) 
a 



^For more general manifolds M, p{x) is smooth on M \ Cut{p). See |19| for precise definitions and statements. 
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Lemma 2.4. {19: . p. 78] Let M he a general n- dimensional complete Riemannian manifold. Sup- 
pose that f : B{{1 + t)R) — > [0, oo) is a non-negative subharmonic function on B{{1 + t)R) (i.e., 
f>0 and Af >0 on B{{1 + t)R) ). Then 

I IWP <4m I (2.18) 

JB(R) ^ JB((1+t)R) 

where C is some universal constant. 

Theorem 2.5 (Gradient Estimate |25[ I19 j ). Let M he an n-dimensional complete Riemannian 
manifold with Ric(Af) > — (n — 1)K, for some constant K > 0. If u is a positive harmonic 
function on M and Br{x) is a geodesic hall in M , then 



IVmI < Cr 



[l+iMy^u onBr_{x), (2.19) 



where Cn is a constant depending only on n. 

Lemma 2.6. ^IQ, p. 15/ Let M be a Riemannian manifold. Then the following holds in normal 
coordinates at x 

A[|VF|2](a;) = 2 ^[a,aji^]2(x) + 2 ^ a,;F(a;)a.(AF)(a;) + 2 Ric(Vi^, VF)(a;). (2.20) 

2.3 Comparison theorem for Jacobi fields and the growth rate of the 
volume form on negatively curved Riemannian manifold 

In this subsection, we only focus on a complete, simply connected, 2-dimensional Riemannian 
manifold M with sectional curvature satisfying —6^ < Km < — < To begin, let O be a se- 
lected point in M, and let expo ■ TqM — > M be the global exponential map at O, whose existence 
is ensured as before by the Cartan-Hadamard theorem. We remark that the tangent space TqM 
is identified with the Euclidean space R^. 



Let (r, 9) be the ordinary polar coordinates on M.^ in the Euclidean sense. That is, the respective 

{x,y) 



induced vectors -§= and -§= are given by 



d_ 

df 

d_ 



i-' + y')- (2.21) 



Then, the geodesic normal polar coordinates (r, 9) on M (as induced by expo ■ ToM — > AI) is 
given by the composite function (r, 9) — (r, 9) o {expo}"^ . 

Let V e ToM = be any unit vector, and consider the geodesic c : [0, oo) -> M with c(0) — O 
and c'(0) — v. Then, we notice that ^^^^^ — c'(r). Next, in order to compute ^^(•^■jj we first 

observe that the ordinary Euclidean -S= is given by 

• -S= — rw, in which w is the unique unit vector in R^ such that the pair {v, w} forms a 

oO rv 

positively oriented orthonormal basis for R^. (Recall that we have the identification ToM = 
R^, so, we may just regard v e ToM to be a vector in R^.) 
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Now, let Y{r) be the Jacobi field along the geodesic c : [0, oo) — ?• M with Y{0) — 0, and 
Y(o) = w. (Recall that Y{r) is a Jacobi field means that the equation 'Vc'Vc'Y + R{Y, c')c' = 
holds along the geodesic ray c{t).) Then, it is well known that [lU] 

Y{r) = {Dexpo)rv{rw) = {Dexpo)rvA ), (2.22) 

On rv 

which implies that ^ — {DexpQ^)c(r)iY{r)). Hence, we have for any / G C°°{M) 

= [{Dexp^^)c(r)iYir))]{f oexpo) = d{f oexpo)rv {{Dexp^'^)c(r){y{r))) 
= {dfUr)iY{r))=Y{r){f). 

This shows that we have the important identity 

d 

— =Y(r). 

dec(r) ^ ' 

With the above preparation, we can now discuss the growth rate of the volume form on a 
complete, simply connected, 2-dim Riemannian manifold with —b^ < Km < — < 0. Under the 
geodesic normal coordinates (r, 9) , the volume form is given by 



d_ 
86 



drdO, 



but sometimes we write | ^ | dpd6 in the case when the distance function p from O is used to 
replace the symbol r. Then, the following comparison theorem in differential geometry is used to 
give us the growth rate of the volume form. 

Theorem 2.7. [Comparison theorem for Jacobi fields \10^ 1 Let M be a simply connected, complete, 
n-dim Riemannian manifold M with —b^ < Km < — < 0. Let c : [0, oo) — > M be a geodesic ray 
in M , and we consider a given Jacobi field Y{r) along the geodesic ray c which is orthogonal to 
c and satisfies |i^(0)| — 1. Then, it follows that we have the following estimate for the growth of 
\Y{r)\ along c, for all r >0. 

1 1 

- sinh(ar) < \Y(r)\ < - sinh(6r). (2.24) 
a b 

Now, by the above comparison theorem together with the identity ^^^^-j — Y{r), we immedi- 
ately deduce that, for a complete, simply connected Riemannian manifold M of dimension 2 with 
—6^ < Km < —a^ < 0, the weight G(r, 6*) = | ^ | of the volume form | ^ | drdO on M is at most of 
the order ^ sinh(6r) . This is used in Section 4 to estimate the integral of a certain non- negative 
function on M. Finally, we also remark that, in the special case of a space form H^(— a^) with 
constant sectional curvature — a^, we have exactly \-^\ = |^('')| = isinh(ar). 



2.4 Geodesic balls, cones and the geometric boundary S{oo) 

In this subsection, we will consider only simply connected, completed, n-dimensional Riemannian 
manifold M with sectional curvature satisfying —b^ < Km < > 0. Recall that the classical 
Cartan-Hadamard theorem ensures that the geodesic normal coordinate at any point O in such a 
negatively curved manifold M is globally defined, and hence M is diffeomorphic to M". 
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Using the distance function from the Section I2.2[ define a geodesic ball in M with radius R 
and centered at x by 

Br{x) = {y G M : p,{y) < R}. 
Next, let O e M and v £ TqM. Define the cone about v with angle 9 by 

Coiv,9) ^ {y e M : Z{v,0^) < 9}. 

Finally the geometric boundary, the sphere at infinity 5(oo) is 

5(00) — the set of all geodesic rays from O, 

which can be canonically identified with the unit sphere in TqM: 6*0(1) = {v G TqM : \v\ — 1}, 
so that every unit vector v € TqM can then be regarded as a point in S{oo), and we can simply 
write V 6 S{oo) (See [HI p.36], [HE]). 

2.5 Bounded harmonic functions on negatively curved manifolds 

Anderson [T] and Sullivan |21| independently, and using different methods, proved the following 
theorem. 

Theorem 2.8. 121}/ Let M be a simply connected, n-dimensional, complete Riemannian man- 
ifold with sectional curvature Km satisfying —b^ < Km < —a^ < 0. Then there exists a unique 
solution u e C°°(M) n C"(A/) to the Dirichlet problem 

Aw = in M, 

A simpler proof is also presented in the comprehensive work of Anderson and Schoen [2] , and 
it is also exposed in [19j . The main idea there is to construct two barrier functions and use the 
Perron's metho4f|. This in turn is accomplished in three steps, whose conclusions we use in the 
proof of Proposition 13. ip . which is a crucial tool for our result. Therefore, we give a brief outline 
of the proof in [191 P- 37], and list the needed conclusions: 

Step 1) Extend the function (/> to all of M and show 

sup \<jy{y) - 0(x)| < Co WW^^^^s^oo)) ^"""^"^ (2-25) 
yeBx(i) 

To extend to all of M we pick a base point O G Af and use the the geodesic normal polar 
coordinates (r, 0) at O to define 

0(r,6') = 0(6*), for all r > 0. 
Lemma [2.31 is then used to show (|2.25p . 
Step 2) The Laplacian of the average of (j) has an exponential decay. More precisely, let 

^^^^ ^ ! x{pl{y))<t){y)dy 
S x{pl{y))dy 

where x is a- standard cut off function. Then it can be showed 

|A0(x)| <Co||0'||^^(5(^))e-"''(-), (2.26) 
where p{x) is the distance function defined in Section [2.21 

•^We note that the proof in [l] also reUes on the Perron's method. 
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Step 3) Show there exists a > and 6 smaU enough such that 

A[0 - ae-^P] > and + ae-^P] < 0. (2.27) 

Then by Perron's0 method there exists a harmonic function F such that 

Ti^-ae-^P <F Kli^ + ae-^P. (2.28) 
The boundary conditions are easily checked. 

2.5.1 Constants a and 6. 

Constants a and d from Step 3 play a very important role in our proofs. Therefore we take 
some time now to discuss a and S and and how they relate to the function (f> and the curvature 
—b^ < Km < — a^. We emphasize that this exposition is completely based on (THJ p. 40] although 
the details of (|2.30p and (|2.31l) below were not exposed there. 

First start with some i5 > to be specified later. Using (|2.13|) we then observe 

A(e-'5''(^)) = Se-^P^'^HS - Ap{x)). (2.29) 

Also, by (|2T^ 

Ap > {n ~ l)a ■ coth(ap) > {n — l)a 

Next, one has to choose sufficiently small S and sufficiently large a so that (|2.27p does indeed hold. 
Let 5 < a, then the first equation (|2.27p is obtained as follows. By (I2.26P and (|2.29p 

A[^ - ae-^P] = A^ - a6e-^P{S - Ap) 

> "Coll^'lUe""" + a5[{n - l)a - 6]e-^P (2.30) 

> {aS[{n ^ l)a - 6] - CoW\\^}e-'P, 

where S < a is used to obtain the last line. Similarly 

A[^ + ae~^P]<CoU'\\^e-''P -aS[in-l)a-d]e-^P 

< {CqU'Woo - oc6\{n - l)a - SWe-'P. ^^''^^^ 

So, for any 5 < a, we choose a = ^Yn-Vja-^] ■ Note, in order to guarantee a > for n = 2, we 
need b <a and not just b <a. Then (j2.27p follows as needed. 

In addition, besides b not being too large, we eventually need b not to be too small. More 
precisely, when we want to obtain that Vi^ is in L^(H^(— a^)), we impose additional condition 
I < 5. Then by the discussion in the Section 12.31 the exponential decay of VF obtained in 
Proposition 13. II below will be sufficient to give || VFj|i2(H2(_a2)) < oo. 

Similarly, when we want to obtain that Vi^ is in l?(M\ where M is complete, simply connected 
2-dim manifold with sectional curvature satisfying —\? < Km < —a?, where a, 6 > and | < a, 
then there we require (5 > |. 

2.6 Hodge Star operator and Hodge Laplacian 

Let d denote the exterior differentiation operator, which sends k forms to fc + 1 forms. As is 
well-known, d satisfies 

ddijj = for any k — form w. (2.32) 

^See [l] for the application in this context or [6]. 
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Its dual operator, d* , is given by 

d* = {-!)'' ***d*, (2.33) 

where * is the Hodge * operator and k comes from d* acting on some given fc-form (see for example 
[18]). We note d* sends k forms to A; — 1 forms. However, the only main fact that we need to know 
in this paper, besides ()2.32p and p.35p below is the definition of the Hodge Laplacian: 

-Aoj = {dd* + d*d)oj. 

When A acts on a function F, then the expression simplifies to 

-AF^d*dF. (2.34) 

So for example if we have a function F that is harmonic, and if we define a 1-form U by 

U^dF, 

then it is very easy to see that J7 is a harmonic 1-form since 

{dd* +d*d)U ^ dd*dF + d*ddF 

= dd*dF by 
= d{d*dF) 
= 0, 

where in the last line we used the fact that F is harmonic and (I2.34|) . The construction of our 
non-unique solution relies on this simple observation. 



3 Exponential decay of the gradient of a bounded harmonic 
functions on negatively curved manifold 

The main result of this section is the following proposition. 

Proposition 3.1. Let M be an n-dimensional complete, simply connected Riemannian manifold 
with sectional curvature satisfying —b^ < Km < —a^. Let 4> G C^{S{oo)) be any boundary data, 
and F £ C°°{M) n C'^(M) be the unique bounded harmonic function on M with F\s(oo) — 4'- Let 
6 < a. Then, the following inequality holds 

\VF\{x)<Co{l + -.. \. ^AU'IU-'"^^^ VxSM, (3.1) 

o[[n — l)a — d\ 

where Cq depends only on a, b, and n. 

Remark 3.2. In the proof, we use Cq to denote a generic constant, which may change from line 
to line, but it always depends only on a, b, and the dimension n of the Riemannian manifold. 

Proof. Given if) G C^{S{oo)) by Theorem 12.81 there exists a unique harmonic function F G 
C°°(A/) n C°(M) with F\s(oo) ^ (t>- By (12:281) we also have 

'^-ae~^P < F <'$ + ae-^P, (3.2) 

where 4> is as in (|2.26p and S < a and a — g^l^^f-^^^^^^ as discussed in Section [2.5.11 Let x G M 
and consider two cases. 
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Case 1: p{x) > 1. 

Consider a ball B^{1). By ^ 



osc F := sup F — mi F 

Next since infB^(i)[(/) — ae~^''] > inf^^ji) (j) — a ■ sup^^^^^ e^^'' and sup5^(i) e~^P = e~^^''^^^~^\ it 
follows that 

osc F< osc + 2ae"'^(''(^)-i\ 

B,(l) 

which implies that the following inequality is valid on Bx{l) 

0<F- M F < osc F < osc + 2ae"'^(''(^)-^^ 

B^{1) B,(l) 

Now, notice that since F — infB^(i)F is a positive harmonic function on Bx{\-) we can apply the 
gradient estimate, Theorem 12. 5[ to deduce the following inequality for any y S Bj.{^) 

\VF\{y) = \V[F{y) - infB^^ijF]\ < Co{l + a)[F{y) - 

< Co{oscB,(i)0 + 2ae~*(''(")-i)} 
<Co[CoU'\\oo + 2ae']e-'P^-^ (3-3) 

So, in particular, if we choose y = x in the above inequality, we have the important conclusion 

\^F\{x) < Co{CoU'\\oo + er/^°"5"°° ^. e^}e-^><-) VxeM-S„(l). (3.4) 

o[(n — l)a — oj 

We now consider the case of x E Bo{l). 
Case 2: < 1. 

Here we have < e^''''*^^^ and sup^^^j^) e^*'' = 1. Hence 

osc F < osc (j) + 2a sup e^*'' 

< Co||(^'||ooe-'^'' + 2a 

= Co||<^'||ooe-'"' + 2ae"e-" 

< [Co||0'||oo + 2ae^e-'^'' 

< [Co||(/.'|U + 2ae^e-*'' 

o[[n — Ija — oJ 

Next, as in Case 1 we can apply the gradient estimate. Theorem 12. 5[ to — iufg^f^i-^F to obtain 
for any y G B^i^) 

\VF\{y) - \V[F{y) - inf F]| < Co(l + a)[F(y) - inf F] 

< Co osc F , > 

Bx(i) (3.5) 
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By taking y = a; in the above inequality, we deduce 

NF\{x) < Co{Co\W\U + .r/^°'n "°° ,, e°}e-^^^-) e Bo{l). (3.6) 

o[[n — l)a — d\ 

By combining p.4p and p.6p we have that p.ip holds for all x e M as needed. □ 

By p.ip and the discussion in Section 12. 3[ we immediately have the following corollaries 
Corollary 3.3. In addition ifS>^, then 

Corollary 3.4. Let M be a complete, simply connected 2-dim manifold with sectional curvature 
satisfying —6^ < Km < —a"^, where a,b > and | < a, if ^ < S < a, then 

1 2 1 1 • • 

4 The proof that ||V|VF| ||Li(H2(-a2)) is finite 

In the proof of ||V|Vi^|^ \\L^(M^(-a^)) < oo, we need the assistance of the following geometric 
lemma, which is itself a consequence of lemma 12.31 

Lemma 4.1. [Covering Lemma] Consider M to be a simply connected, complete 2-dimensional 
Riemannian manifold with sectional curvature —b^ < Km < —a^ . Let O be a selected base point 
in M , and let p be the distance function from O. Then, there exists some sufficiently large uni- 
versal constant Rq > 2 such that the following assertion holds 

For any given R > Rq, if we take the positive integer N{R) — [ ] + 1 = [ttg^^] + 1 (here, 
the symbol [A] means the largest integer N E Z with N < X), then, we can pick a list of vectors 
vi,V2,V3, ....wjv(fl) G '5'(oo), which are evenly distributed on the circle ^(oo) in such a way that we 
have the following inclusion 

{xeM:R-l< p{x) < i? + 1} C Ut^f)s,„^(fl) (3(1 + \)), (4.1) 

where for each 1 < N{R), Cy. : [0, oo) — > M is the geodesic ray of unit speed with Ct,;(0) = O, and 
<(0) = z;,. 

Remark 4.2. In words, what the conclusion of the lemma is saying is that if we consider an 
annulus in M with inner radius R — 1 and outer radius R+l, where R is big enough, then we can 
cover it by N{R) geodesic balls centered at (i?)with radius (3(1 + -j)). 

Proof. To begin, let us select a base point O in M, and let p be the distance function from O. By 
Lemma 12.31 there exists a sufficiently large universal constant Rq > I such that for any two points 
xi, X2 in M with R = \Oxi \ = \Ox2 \ satisfying R > Rq, we have the following 

2R+-{log9-l) < \xTx^\ < 2R+^{\og9 + l), (4.2) 
a 



where 9 = Z{Oxi,Ox2)- 

From now on, we use the universal constant Rq 
V S S'(oo). We then consider the geodesic ray 



— Rq + I. Now, choose any R > Rq, and let 



Cy : [0, oo) M, 

c„(0) = O, and Cy{0) = v. 
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Now, we consider the universal angle 9^-^^ = e and the sector To{v,6'^^^;R— l,i?+l) defined 

by 

To{v,9^"'>;R-l,R+l) = {x e Co{v,9'^"''>) : R - 1 < p{x) < i?+ 1}, 

where the cone Co{v,9^^^^ was defined in Section [^Hl Our goal is to prove that To{v,9^^^;R~ 
1) C Bc„(H)(3(l + i)). To this end, let x £ To(w, 6l(^); i? - 1, i? + 1). Then, p{x) = R + X, 
for some A £ [—1,1]. By the triangle inequality, we have 

\cyiR)x\ < \cy{R)cy{R + X) \ + \c^{R + X)x\ < \X\ + \c„{R + X)x\. (4.3) 

But from (14. 2p with xi — Cy{R + A) and X2 = x, it follows 

\c^{R + X)x\ < 2(i?+ |A|) + |{log[Z(Oc,(i? + A),Oi)] + 1} 
< 2i? + 2+|[log((?(^)) + l] 



= 2i? + 2 + |{log(e-''^) + l} (4.4) 



= 2R + 2 + y{-bR + l) 
b 

Hence 

\c^{R)x\ < \X\ + \c,{R + X)x\ < |A| + 2(1 + i) < 3(1 + i). (4.5) 

This shows that every x eTo{v,9'-^^;R-l,R+l) must lie in the geodesic ball B^^^r) (3(1 + i)). 
To conclude the proof, we just take the integer N{R) = [ ^^-hR ] + 1- Then, we can select some 
evenly distributed vectors wi,U2, ...,ujv(fl) S ^(oo) such that 

{a; e M : i? - 1 < p{x) <R+l} = U^J^ho{v„9'^^'^; R - 1, i? + 1). (4.6) 

bmce we already know that, for each 1 < i < N{R), we have Toiv.i,9'-^^- R - l,i? + 1) C 
i?c„.(fl)(3(l + ^)), in which Cy. : [0,oo) — > M is the geodesic with c«.(0) = O and c^.(O) — Vi, it 
follows at once from relation (|4.6I) that 



{xeM : R-l< p{x) <R+l] = U^J^''To{v^,9^^-'>;R-l,R+l) 

cufjf)s.„^(,)(3(l + i)), ^'-'^ 

as desired. 

□ 

With the help of the Covering Lemma 14.11 and Lemma 12.41 we can now prove the following 
fact. 

Proposition 4.3. Let a,b > satisfy ^b < a < b, and let M be a simply connected, complete 
2-dimensional Riemannian manifold with sectional curvature —b^ < Km < —a^. Then, for any 
bounded harmonic function F € C°°{M) f] C"(M), which arises from - boundary data (j), it 
follows that /jjV|Vi^p| < oo. 
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Remark 4.4. We note that in the proof of Proposition 14.31 it is not necessary for us to obtain a 
uniform bound of Jj^ |V|VFp| < oo in terms of, say, H^'Hoo- AH we need is just to confirm that the 
integral /^^ |V|Vi^p| is finite, because this is aheady enough to ensure that /^^ dw{V|VFp} — 0. 



Proof. As usual, we begin with a bounded harmonic function F G C°°{M) n C°(M) such that 
^ls(oo) ~ 'f'' sufficiently large universal constant as determined in Lemma |4. II 

Since F is smooth on M, in order to prove that /^^ |V|VFp| < oo, it is sufficient to see that 
Im~Bo{R) I^I^^I^I ^ ^'-'^ some large R > Rq, where O is a selected base point in M. 
We first write 

|V|Vi^n < |V{|Vi^p + Ae-^^P}\ + A\ye-^^P\ 

= |V{|VF|2 + Ae-^^P}\ + A{2S)e-^^P. ^^'^^ 

We estimate the first term on the right. First, by p.20p 

A[\VF\^]{x) = 2[d^djF]^{x) +2d^F{x)d^{AF){x) +2Ric(yF,\7F){x). 

Since AF = it follows from the above formula that 

A[|VFp] > 2Ric(VF,Vi^) > -26^1 VFp. (4.9) 

To proceed further, we take S to be any fixed choice of positive number within the range 
I < 5 < a (i.e., we choose such a S once and for all) and by Proposition 13. II we have 

\VF\ <CaMW Woce-'", (4.10) 

where the constant Ca.b depends only on n = 2 and a, and b. Notice that our fixed choice of 
6 G (|,a) automatically satisfies the condition S > ^, due to the fact that b > a. Next, since 
(5 > |, we can choose some sufficiently large positive number depending on 5, R{6) > 2, such that 

Next, we have to find some A > large enough, and some sufficiently large radius R such that 
the function |Vi^p + Ae~^^P will be subharmonic on — Bo{R)- To achieve this we use (|2.15p 
(with the condition Km > ~b^) and observe 

ApX{^>«m} < b ■ coth(6p)X{p>flffij < (& + -)X{p>!u^} < &(1 + -^)X{p>3my (4.12) 

Hence from (gl]), (|4T0l) and ([2:29| 

A{\\/Ff + Ae-2^''}X{^>K(^j > {-2b'Cl,U'\\le~''P + Ai2S)e-''''[2S - Ap]}x{^>«^j 

1 . ' (4.13) 



> 2{AS[2S 6(1 + — -)] - b'Cl,U'\\lo}e'''%p^iM,^ 



Notice that we definitively have 26 — 6(1 + - j^) > 0, thanks to our choice of R{6) which ensures 
the survival of the second inequality sign in 14.111 Next, we just take 



<5[2<5 - 6(1 + ^)] 
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With this choice of A, it follows from (|4.13p that 



A{|VF|2 + Ae-^'p}x,„:,n^, > 0. (4.15) 



4p> 

That is, the function |Vi^|^ + Aer'^^f is subharmonic on M — Bo(p^\ So, we may apply 
Lemma to \V ^ Ae-'^^P and deduce that for any geodesic baU B^(6(l + -i)) C M-Bo{^), 
we have 



/ |V{|VF 




|V{|VF|2+^e-2A>||2 I 
(3(1+^)) J (4.16) 

1 

\VF\^ + Ae-^^Pf 

/Bx(6(l + i)) J 

<C(a,6;||0|U)e-2^''(^), 

where we again used (I4.10p to go to the last line. We further remark that, in the above 
estimation, we have implicitly employed the volume comparison theorem in differential geometry 
which says that Km > —b^ implies that the volume of any geodesic ball Bx{6{l + -j)) in M is 
bounded above by a universal constant C'b (such a universal constant Cb which serves as the upper 
boimd is indeed the constant volume of any geodesic ball with radius 6(1 + ■^) in the space form 

H2(-62)). 

Now, let us take Kq to be a sufficiently large positive integer for which Kq > max{Ro, ^^^^ + 
6(1 + -j)}, where Rq is the sufficiently large universal constant determined in Lemma [4.11 Then, 
by Lemma [4.11 for any positive integer k > Kq, if we take the positive integer N{k) — [ire''''] + 1 
then, we can pick a list of vectors Vk,i, Wfc,2, "^fc.s, ■■■■Vk,N{k) G '5'(oo) in such a way that we have the 
following inclusion 

{xeM:k-l< p{x) < fc + 1} C ut?^S,„^ jfc)(3(l + i)). 
By combining inequality (j4.16p and the above inclusion, it follows that 

/ |V{|VF|2 + ^e-2^P}| < V / |V{|VF|2 + Ae^25p|| 

J{p{x)>Ko} k=Kohk-l<p{x)<k+l} 

oo N(k) 

< E E / \V{\VF\^+Ae-^'P}\ 



k=Ka i=l + 
oo N(k) (4.17) 

< E E^(«'^'ii<^ii-)^"''' 

oo 

< E C{aA\moo){[iTe'"'] + l}e-^"' 

k=Ko 

< CO. 

We note that to obtain the last inequality X]fe!L_fCo C!{a, h, ||0||oo){[7''e'''''] + l}e^^*'^ < oo, we use the 
fact that our fixed choice of 5 lies within the range ^6 < J < a , so that [25 — 6) > is automatic, 
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and hence ZZko e''''e~^^'' = ^Zko e-'^^"^)^' < oo. 

Next, we notice that the volume form on M (with respect to the geodesic normal polar co- 
ordinates {r,6) about O) is in the form of G{r,d)drd6, where G{r,9) < i sinh(6r) thanks to the 
comparison theorem 12.71 for Jacobi fields (with Km > —b'^) ■ Hence, it follows again from 26 > b 
that 



g-25p^ / / e-^^'^G{r,9)drd9 

M Jo Jo (-^-^g^ 



</ / e-'^^''-smh(br)drd9 < oo. 
Jo Jo b 

Hence by g^, (|iTfl) and 

/ \y\yFf\< [ \\7{\\7Ff +Ae-^^P}\+A{26) [ er^^P < oo. (4.19) 

J{pix)>Ko} J{pix)>Ko} J{p(x)>Ko} 

Since [VlVF^j is continuous in M, bv (|4.19p we must have /^^ |V|Vi^|^| < oo as needed. □ 

5 Finite Dissipation 

We begin with two propositions, which help us establish the energy inequality (|1.7p . First, using 
g(-, •) on T*M ® T*M defined in Section [21 we can consider for each 1-form 9, two non-negative 
valued functions g{Def9,Def9) G C°°{M), and g{V9,y9) € C^{M). We have the following 
relationship between them. 

Lemma 5.1. For any given n-dimensional complete Riemannian manifold M , we have 

g{Def9,Def9)<g{V9,V9). (5.1) 

Proof. Let p 6 M, and consider the geodesic normal coordinates (x^, a;^, x") about the point 
p, so that the natural frame di\p,d2\p, ...,dn\p (induced by the geodesic normal coordinates) at 
the point p is orthonormal, and that the Christoffel symbols F^-^, (induced by the geodesic normal 
coordinate) vanish at the point p. Hence, for any 1-form 9 — 9jdx^ we have 9j.kip) = dk9j{p). 
So, by (I2H) it foUows 



g{Def9,Def9)\p = ^giY^{d,9j + dj9,)ip)dx' ® dx^p,J2idkOi + di9k){p)dx'' ® dx'\p) 

i,j k,l 

= i E + dj9,){p){dk9i + di9k){p)5'H^' 

i,j k,l 

= \ E(^^^^ + dA)iPm9, + d,9i){p) 

id 

= Y,[d,9,{p)f. 



On the other hand, by (|2.12p the Christoffel symbols F' j. vanish at p, so it follows that 9\p 
J2i,j dt9j{p)dx'- (g) dx^p. Hence 
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hi 

and 

Since p € M is arbitrary in the above argument, it follows that the above inequality is valid for 
all points in M as needed. 

□ 

Proposition 5.2. Let a,b > and and such that < a < b. Let M be a simply connected, 
complete 2-dimensional Riemannian manifold with sectional curvature ~b^ < Km < —a^. Let 
4> G C^{S{oo)) be any given boundary data, and let F S C°°(A-f) n C°(M) be the unique bounded 
harmonic function on M with -F|s(oo) — 4'- Then the following holds 

[ g{V{dF),^{dF)) = - [ Ric{VF,\/F)<b^ [ |VFp. (5.2) 
Jm Jm Jm 

Proof. Let F G C°°(il/)nC°(M) be the unique bounded harmonic function on M with i^|s(oo) = 4'i 
where (f) G C^{S (oo)) is some given boundary data. 
Then again by (p:20)) 

A(|VF|2)(a;) = 2 ^ [^^^JFf{x)+2 ^ ^^F{x)^^{AF){x) + 2mc{\/F,VF)ix), (5.3) 

1<«J<2 l<ij<2 



where di , 82 is the natural coordinate frame induced by the geodesic normal coordinates about the 

- jk 



point X. So again by (|2.12p Christoffel symbols F' j, vanish at x, and we have 



J2 mFnx)^giV{dF),VidF))U. 

l<ij"<2 

Also, as before we use AF = in (|5.3I) to obtain 

^A{\\7F\^){x) ^g{V{dF),V{dF)){x)+Ric{\7F,VF){x) Vx G M. 

Now, for each positive integer k > 1, consider a smooth function ipk ■ [0,oo) — )• IR, which 
satisfies X[o,2'^] < 4'k < X[o,2'=+i]j and IV'^I < Now, let O be a selected base point in M, and 
let p be the distance function from O. Then, by multiplying the above equality by the cut off 
function ipkip^) (which is compactly supported in Bo{3k)) and integrating over M, we yield the 
following equality 

/ ]-div{V\VF\^)Mp')^ I MP^)9i^{dF),V{dF))+ [ Vfe(p2)Ric(VF,Vi^). (5.4) 
Jm ^ Jm Jm 

But since ipkip'^) is compactly supported in Bo{'iR), it is plain to see that, for every fc > 1, we 
have 

1/ \d^v{V\VF\^)MP^)\ = \\ I 2V^l.(p2)Vp.V(|V^^p)| 
Jm 2 2 Jm 

< ttt: / |V(|V^^|2)|. 



2k 



m 
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Since, according to Proposition 14. 3[ we have /^,^ |V(|Vi^|'^)| < oo, it follows from the above 
inequality that 

lim / idw(V|VFnV'fe(p2) = 0. (5.6) 
On the other hand, by the monotone convergence theorem, we have 

hm / MP^)9{^{dF),V{dF))^ [ g{V(dF),V{dF)), (5.7) 

and that 

lim / i^fe(p2).(_Ric(VF,Vi^))= / (-Ric(Vi^,VF)) (5.8) 
As a result, by taking the limit of each side in equality (|5.4p . we get 

0== / g{V{dF),V{dFj)+ [ Ric(VF,VF). (5.9) 
Jm Jm 

That is, we have 

/ g{V{dF),V{dF)) = - [ Ric{\/F,VF) < [ \\/F\^, 
Jm Jm Jm 

in which the last inequality follows from the fact ^^JU\i that Ric(VF, Vi^) = /CmIVF^ > 
-62|VF|2. □ 

Corollary 5.3. Let a > and let (f> £ C^(5(oo)) be any given boundary data, and let F g 
C°°(]HI^(— a^))nC"(H2(— a^)) be the unique bounded harmonic function onIP{—a'^) with F\s(oo) = 
(j). Then the following holds 

I g{^{dF),V{dFj)^a^ [ \\/Ff . (5.10) 

Proof. This is immediate from equation (|5.2p since Ric(Vi^, VF) = -a2 \VF\^ . □ 

6 Proofs of the main results 

First we establish the following lemma. A simpler computation in normal coordinates could also 
be done in the same spirit as the computation in the Euclidean space. However, we present a 
different proof below due to its intrinsic nature. 

Lemma 6.1. The following identity is valid for any smooth function f on any given n-dim Rie- 
mannian manifold M . 

Vv/d/ = id|d/|'. (6.1) 

Proof. First, for any smooth vector field AT on a Riemannian manifold M , and for any smooth 
function / on Af, we write 

^if) = {X,d,f)xM^T*M- 

Next, recall that the Lie bracket [AT, Y] between two vector fields X and Y, is itself another vector 
field, and is characterized by [X,Y](f) = X{Y{f)) - Y{X{f)), for any / G C°°{M). 
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Now, for a given smooth function / on a Riemannian manifold M , we consider the gradient 
field W — V/, which means that W* — df. Then, by (j2.ip it is plain to see that 



W{f) = {W,df)TAmT'M = 9{W,W) = \W\^ (6.2) 

Next, we have the following identity for any smooth vector field X on M, due to the fact that the 
Levi-Civita connection V on TM is compatible with the Riemannian metric g{-, ■) on M we have 

^X{\W\^) - ^{X,d{\W\^))TM^T'M = g{^xW,W), (6.3) 

and 



{X, \7wW*)tm»t'M - {X, {VwW)*)TAmT'M by jM]) 
^g{X,VwW) by dSD) 

= W{g{X, W)) - giVwX, W) by compatibility 
= W{{X, df)TKmT'M) - g{VwX, W) 

^W{X(f))-g{VwX,W). (6.4) 

But due to the torsion free property of the Levi-Civita connection V on TM, which says that 
-Vm/X = [X, W] - WxW, we have 

-g{^wX,W)^g{[X, W],W) - gi^ xW,W). 

Hence by ([O)) 

{X,VwW*)TivmT'M = W{X{f)) + g{[X, WIW)- g(VxW, W) 

= W{X{f)) + [X,W]{f) - ix(|M^p) by ^ and ^ 

= W{X{f)) + X(Ty(/)) - W{X{f)) - \x{\w\^) 

where the last equality follows since X{W{f)) = X{\W\^) by (|6.2p . In conclusion, by using (|6.3p 
again, the following equality holds for any smooth vector field X on M 

< X.ywW* - ^di\W\^) >TAmT'M= 0, 

which means the same as saying that Vv/d/ - |d(M/P) = as needed. □ 

6.1 Proof of Theorem [US 

First we show existence and the lack of uniqueness. 

Existence and Non-uniqueness: For convenience we recall the Navier-Stokes equation on 
H2(-a2). 

dtU* - AU* +VuU* -2Ric{U*) +dP = 0. 

d*U*=0 (N-Sh^(-=)) 
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Now, let (f) G C''\S{oo)), then by Theorem 12.81 there exists a (unique) harmonic function 
define a solution {U* ,P) to be 



F e C°°(H2(-a2)) n C"(H2(-a2)) satisfying F\g^^^ = 0. We let our initial data uq = dF, and 



U* = P = -dtip{t)F - i + 2a^F (6.5) 

where is any function satisfying 

^'^{t) + 2a^ [ ^p^{s)ds <i^^{0). (6.6) 

For example, we could let ilj{t) = exp(— ^) for some A > 2a^. 

First we show {U*,P) solves N-Su2(_a2). This is very easy by the preparations we have done 
in Section!!^ Indeed, by AU* = 0, and we observe that by Lemma [Ol and d^JT]) . 

dtU* +VuU* + 2Ric(C/*) = -dP. 

It is also very easy to see that U* is divergence free since by definition of U* and (|2.34l) this is 
equivalent to F being harmonic. 

Global energy inequality: Recall we want to show 



[ \U*{t,x)\^ + 2 [ I g(DeiU*,DeiU*){s,x)ds< [ |uo|^ 

JH2(_a2) Jo JH2(-o2) JH2(-a2) 



Thanks to Lemma [5. II we have 



/ \U*{t,x)f + 2 f f g{DefU*,De{U*){s,x)ds 

JH2(-a2) Jo JH2(-a2) 

< / \U*{t,x)f + 2 f f g{VU*,WU*){s,x)ds 

JH2(-a2) Jo JH2(-a2) 



By Corollary 



/ g{V{dF),V{dF))^o? f \VF{t,x)f 

JH2(-a2) JH2(-a2) 



SO 



g{VU*,VU*) = t/j^it) \VF{t, x)\ 

H2(-a2) JH2(-a2) 



Hence by and (|6j 



ft 

2 



\U*{t,x)f + 2 / / g{VU*,yU*)ds 

H2(-a2) Jo Jh2(-o2) 

^Jj\t)\VF{x)\^ + 2a^ f f V'(s)|VP(x)|^ds 

H2(-a2) Jo JH2(-a2) 

iP\t)+2a^ f\/{s)ds] f |VP(x)|' 

Jo / JH2(-q2) 

< tp^iO) [ |VF(a;)l^ 

JH2(-q2) 

\uo\ , 



H2(-a2) 



as needed. 



26 



6.2 Proof of Theorem 11.61 

The proof is very similar. Therefore we just give a brief sketch. 

Existence and Non-uniqueness: Again, for convenience, we recall the modified Navier-Stokes 
equation on M. 

dtU* - AU* +yuU* + dP ^ 0, 

,*TT* n (N-Sm) 

d*U =0. 

Let uo — dF, and 

U* = i:{t)dF P = -dtijj(t)F - ^ \dF\^ (6.7) 

with ijj{t) — exp(— ^) for some A > 26^.. Then as before we can see the equation is satisfied. 
Global energy inequality: Bv l5.ll we need to establish 

\U*{t,x)\^ + 2 I I -g(yU\VU*)ds< I \uq\^ . 

M Jo Jm Jm 



By Proposition 15.2 



so 



Hence 



/ 5(V(dF),V(di^)) <62 / |VF^ 
Jm Jm 

f g{VU*,WU*)<b^ f i^^{t)\\7F{t,x)\ 

JM JM 



\U*{t,x)\^ + 2 [ [ g{VU*,VU*)ds 
M Jo Jm 

<f ■tl;^{t)\VF{x)f +2b^ f f \ij{s)\/F{x)fds 
Jm Jo Jm 

i^^{t)+2b^ f V^(s)ds) f |VF(a;)|^ 
Jo ) Jm 

<V2(o) f |VF(a;)|' 
Jm 



M 



as needed. 



6.3 Proof of Corollary [Tl and [TT] 

In the recent paper Koch, Nadirashvili, Seregin and Sverak studied Liouville thorems and 
their consequences for the Navier-Stokes equations. One of the results is 

Theorem, fllf Let u be a bounded weak solution of the Navier-Stokes equations on x (— oo,0). 
Then u{x,t) = b{t) for a suitable bounded measurable function b : (oo,0) — > R^. 
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The three dimensional problem is more difBcult. Nevertheless Koch, Nadirashvili, Seregin and 
Sverak are able to obtain corresponding results for the axi-symmetric equations with no swirl. 
What Corollaries 11.41 and 11.71 show that in the hyperbolic setting we can have bounded solutions 
in both two and three dimensions (in fact, for any n > 2) that are not functions of time only. The 
nontrivial bounded solutions we choose are in the form of (|6.5p and (|6.7p for N-S/fn(_^2) and (|1.6p 
respectively, where we drop the condition (j6.6l) . which is only needed to show the global energy 
inequality. It would be interesting to find out whether or not these are the only bounded solutions 
of N-Sjj„(_a2) and (HIS]) . 

Here we also mention the result of Galdi 5 , which states 

Theorem, For the steady Navier-Stokes equation on M.^ whenever the solution satisfies the 

9/2 

finite dissipation property and u S LJ , then u must he a trivial solution, i.e. u is constant. 

9 /2 

We note that in our case, we have a nontrivial solution, which belongs to LJ , but at this time 
we cannot say whether or not there exist nontrivial solutions in three dimensions that also satisfy 
the finite dissipation property. 

Proof of Corollary II .41 and 1 1 . 71 Let -0 be bounded in (|6.5p and (|6.7p and without the condition 
(|6.6p . It is obvious from the proof of Theorems 11.21 and 11.61 that (|6.5p and (|6.7|) are the solutions 
of N-S/fn(_a2-) and (|1.6p respectively since this is independent of the dimension of the underlying 
manifold. Hence, we only need to verify that U* is L°° bounded. That can be checked in more 
than one way as follows (note we also do not need the exponential decay of the gradient of F). 
Since F e C°°(M) H C"(Af), \dF\ < cxd or we could use the much more sophisticated tool of the 
gradient estimate. Theorem 12.51 
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